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TTTT scattering amplitudes constrained by Roy's equations*^ 
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The scalar-isoscalar, scalar-isotensor and vector-isovector tttt amplitudes have been fitted simultaneously to 
experimental data and to to Roy's equations. Resulting tttt phase shifts up to 1600 MeV and near threshold 
observables have been analyzed. Only the amplitudes fitted to the "down-flat" set of phase shifts in scalar- 
isoscalar wave fulfill crossing symmetry conditions and can be regarded as physical. 



1. Introduction 

After the ne-w analysis of the tt" — > 7r+ tt^ n 
reaction on polarized target in 1997, the long 
standing t-wo-fold "up-do-wn" ambiguity in scalar- 
isoscalar TTTT amplitudes belo-w 1000 MeV reap- 
peared as a four- fold one [J . T-wo of these solu- 
tions, called "steep", violate unitarity and can be 
rejected as non-physical |1I2| . The phase shifts 
for the f-wo left "flat" data sets are presented in 
Fig. 1. The biggest differences are in the effective 
f-wo-pion mass range TOtttt bet-ween 800 and 970 
MeV i.e. near the t-wo scalar resonances /o(500) 
and /o(980). The kno-wledge of the phase shifts 
in that region is important for example in the 
determination of the parameters of these scalar 
resonances and it is crucial for -works on disper- 
sion relations for the tttt amplitudes belo-w 1000 
MeV 0. 

2. Analysis of the data 

To determine -whether unitary amplitudes, fit- 
ted to the "flat" solutions, can be treated as phys- 
ical ones we have checked if they fulfill crossing 
symmetry conditions 0]. We have used, as in- 
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put amplitudes in Roy's equations, the scalar- 
isoscalar (SO), scalar-isotensor (S2) and vector- 
isovector (PI) -waves defined in |4I5| . We have 
performed simultaneous fits to experimental data 
and to Roy's equations in the m,r7r range from 
the TTTT threshold to 970 MeV using the test 
for Xtot 
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-where (fil (si) and Aip^ (si) represent the ex- 
perimental phase shifts and their errors, re- 
spectively. The real parts Re fl^{si) — 
77 sm[2(5^(si)]/(2 fci s^^) are expressed as func- 
tions of the TTTT phase shifts S^{si) and have been 
calculated under the assumption that the inelas- 
ticity 77 is equal to 1 in the studied m^7r range. 
The other real parts, denoted by Re /o„t(si) are 
the output values calculated from Roy's equa- 
tions. We took A/ value of 0.5 x 10~^ to obtain 
reasonable fits to Roy's equations. The Nj — 18 
experimental values of the "up-flat" or "do-wn- 
flat" data bet-ween 600 and 950 MeV -were used 
in addition to the six data points taken from [5]. 

The SO--wave amplitudes have been constructed 
by fitting the three coupled-channel model predic- 
tions for phase shifts and inelasticities (see |5]) to 
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Figure 1. Scalar- isoscalar phase shifts [11 for the 
solutions "up-flat" (open circles) and "down-flat" 
(full circles); the round curve indicates the re- 
gion where the "up-down" ambiguity leads to the 
largest differences between the two solutions. 



the data of fTl. In the S2-wave we have used the 
data from solution A of Hoogland et. al. [7] and 
in the Pl-onc the data of Hyams et. al. 

Below 970 McV the following Pade representa- 
tion of the SO phase shifts has been taken: 

where k = ^\/s — Am'^ is the pion momentum 
and aj {j = 1,...,7,9) are constant parame- 
ters. Above 970 MeV and up to 2 GeV our 
coupled channel model amplitude A, fitted to 
the "down- flat" data, and the amplitude C, con- 
strained by the "up-flat" data, were used. These 
two representations, one below 970 MeV and the 
second above were joined smoothly at 970 MeV 
up to their first derivative. In the fits we have 
also used the near-threshold phase shifts calcu- 
lated from the differences — Si obtained in the 



high statistics Kf^4 decay experiment ISl. 

The parameterization of the S2-wave, using 
a rank-two separable potential model, has been 
described in ^ where detailed analysis of the 
present study is presented. 

For the P-wave, from the tttt threshold up to 
970 MeV, we have used an extended Schenk pa- 
rameterization 0: 

tan Si (s) = ^ fc^ 

where A is the P-wave scattering length and Sp is 
equal to the p-mass squared. Above 970 MeV we 
took the if-matrix parameterization of Hyams et 
al. 0- The parameters C and D were chosen to 
join smoothly around 970 MeV the phase shifts 
given by both parameterizations. 

The contributions to Roy's equations from high 
energies (totttt > 2 GeV) and from higher partial 
waves {I > 1), called driving terms, are composed 
of contributions from the /2(1270) and p3(1690) 
resonances and from the Regge amplitudes for the 
Pomeron, p- and /-exchanges. For the /2(1270) 
and p3(1690) we have used the Breit-Wigner 
parameterization with masses, widths and tttt 
branching ratios taken from jJUj. For the Regge 
parts we have used formulae of IP without the u- 
crossed terms. In the driving term the /2(1270) 
resonance dominates in the scalar-isoscalar wave 
and the isotensor and isovector waves are mostly 
influenced by the p3(1690). In the isoscalar wave 
the Regge contributions are more than 10 times 
smaller than the resonance contributions but for 
the isospin 1 and 2 they are of the same order. 

In Fig. 2a and 2b we present results of fits to 
the " down- flat" and " up-flat" phase shifts and to 
Roy's equations (solid lines). In both cases the 
differences | Re - Re f^^ \ are of the order 
of 10^"^ in all three partial waves. The for 18 
points between 600 and 970 MeV was 16.6 in the 
"down-flat" and 46.4 in the "up-flat" cases. We 
see in Fig. 2b that the solid line lies distinctly 
below the "up-flat" data points between 800 and 
970 MeV. In contrary, the corresponding line for 
the "down- flat" case in Fig. 2a is very close to 
experimental data in the same range of TOtttt- In 
order to improve the fit to the "up-flat" data we 
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Figure 2, Solid lines: fits to the scalar-isoscalar 
phase shifts of [1] and to Roy's equations; dashed 
lines: fits to the data points shifted upwards 
(called "upper" and drawn as long-dashed lines) 
and downwards (called "lower" and drawn as 
short-dashed lines) by their errors. Dotted line 
for the "up-flat" solution represents fit to experi- 
mental data only. 



have used constraints given by the good fit to the 
"down- flat" data. Two parameters were fixed by 
choosing the values of the scattering length and 
the slope parameter and two others by the values 
of phase shifts calculated from this fit at 500 and 
550 MeV. A new fit with these constrains gave 
smaller value of xlxp = 13 for 18 "up-flat" data 
points, but simultaneously provided us with an 
enormous value of Xjioy = ^-^ >^ lO''- The phase 
shifts for this amplitude are presented in Fig. 2b 
by the dotted line. It is clear that a simultane- 
ous good fit to the "up-flat" data and to Roy's 



equations is impossible. 
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Figure 3. Real parts of input amplitudes "m" 
(their phase shifts have been presented in Fig. 2) 
and of output ones "ouf calculated from Roy's 
equations; all values are multiplied by 2 A;s~^/^. 



Besides the fits to the "up-fiat" and "down- 
flat" experimental points we have also performed 
flts to data points shifted upwards and down- 
wards by their errors. In these fits the same four 
constraints, just described above, were used be- 
low 600 MeV. Up to 937 MeV in the "down-flat" 
case in Fig. 3, the curves labeled upper "in" and 
lower "in" form a band included inside a band de- 
limited by the lines upper "out" and lower "out". 
All the curves lying inside these bands correspond 
to the amplitudes fulfilling crossing symmetry so 
the " down-flat" data can be accepted as physical 
ones. In the "up-flat" case in Fig. 3 the output 
band lies outside of the input band from 840 to 
970 MeV. It means that in this case crossing sym- 
metry is violated by the amplitudes fitted to the 
"up-flat" data. 
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3. Fit to "down-flat" solution with uniform 
parameterizations 

Having chosen the "down-flat" solution as the 
physical one we have performed another fit with 
different but uniform parameterizations in SO-, 
S2- and Pl-waves separately (we have called this 
fit uniform fit). We have used our three cou- 
pled channel model for the SO- and S2-waves 
and constructed an analogous one for the Pl- 
wave. The uniform fit has been done in wider 
m,r7r range, than that described in section 2 (here- 
after called nonuniform fit) . In the SO- wave case 
the fit was done from the tttt threshold to 1600 
MeV and to 1250 MeV for the S2- and Pl-wave 
cases. In table 1 we present values of func- 
tions defined in equations |^ and ^ for both 
fits. The values of ^-re bigger in the uniform 
fit due to a larger TOtttt range where the fit has 
been performed and to a smaller fiexibility of the 
model. 



Table 1 

Values of functions for nonuniform fit (first 
line) and uniform fit (second line); in parenthe- 
ses is number of experimental points; for Xtot 

the 

second number is the total number of free param- 
eters used in the fit. 

Xioi \exp \Roy 

43.9 (80/13) 37.4 (44) 6.5 (36) 
247.5 (245/22) 221.7 (179) 25.8 (66) 



In table 2 we compare the values of the near- 
threshold parameters for both nonuniform and 
uniform fits. Significant differences in the Pl- 
wave case are due to the fact that near-threshold 
parameters of this wave are not sufficiently con- 
strained both by Roy's equations and by the ex- 
perimental data. 

4. Conclusions 

The use of Roy's equations allows final elim- 
ination of the long-standing up-down ambiguity 
in the tttt phase shifts in favour of the "down- 
fiat" solution. This agrees with the resent results 



Table 2 

Values of threshold parameters for SO, S2 and PI 
waves for fits to "down-flat" solution. The first 
raw gives the results for the nonuniform fit and 
the second one those for the uniform fit. In the 
SO- and S2-wave cases the ag and 6g are multi- 
plied by and 10^ rn^ respectively and in the 
case of the PI wave by 10^ and 10'^ m^. Pa- 
rameters for the nonuniform fit with error bars 
are presented in 15. 

^ b-'o ai b\ 

0.22 0.25 -3.4 -7. 4. 2.6 
0.20 0.24 -4.4 -8. 5. 7.2 



of the joint tt+tt" and tt^tt" analysis Addi- 
tional constraints from for example Froissart Gri- 
bov sum rules or from ChPT are needed to fix the 
near threshold behaviour of the tttt amplitudes. 
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